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ABSTRACT 


The moments of a region in an image can be used to describe the 
region's location, orientation, and shape. This paper derives the class 
of all possible fonnulas for computing arbitrary moments of a region 
from the region's boundary. The selection of a particular formula 
depends on the choice of an independent parameter. Several choices of 
this parameter are explored for region boundaries approximated by 
polygons. The parameter choice that minimizes computation time for 
boundaries represented by chain code is derived. Finally, two 
algorithms are presented. The first computes arbitrary moments for a 
region from a polygonal approximation of its boundary. The second 
algorithm is optimal for computing low order moments from chain-encoded 


boundaries. 


I. imODUCTIOM 


Many pattern recognition techniques in con^uter vision use 
structural or statistical featiures of regions and their outlines to 
characterize the shape of objects being viewed. An important class of 
statistical features is the set of (p ••■q )-th order moments defined by 


%■ / / x^y*^f{x,y) dxdy (l) 

rnrnm mrnm 

where f(x,y) is a density distribution function. In the context of 
pattern recognition, the moments are computed for a unifrm density 
distribution over a closed region R in the xy-plane. Therefore, f(x,y) 
reduces to 

i l, if (x,y) eR 
0, if (x,y) fR 

In this case, (l) becomes 


m 


pq 


J'J' x^y*^ dxdy 

R 


(P) 


The most common example of the use of moments is to compute the 
centroid (x,y) of a region by 


m 


X = 


m 


10 

00 




y = 


m 


00 


( 3 ) 


where m^^ is the area of the region. Mcanents are also used to compute 
the angle of a region's axis of minimum moment of inertia 0. This quantity 
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defines the region's orientation within a two-fold degeneracy and Is 
given by 


e 



- (“00®20 ” ®lo) “ (“oo“o2 



{k) 


The quantities x, y, and 0 are useful for pattern recognition 
because they specify the position and orientation of regions, defining 
a transformation between Image and model coordinates, and thus allowing 
further analysis of shape features In a standard reference frame. 

A more direct application of moments to pattern recognition Is the 
use of moment invariants to describe objects. These are qxiantities, 
computed in terms of the moments of a region, that are invariant under 
translation, rotation, and scede changes. Hu [l] derives seven moment 
Invariants and Wong and Hall [2] describe an application that uses 
them to recognize objects in aerial scenes. 

If the region is scanned in a raster fashion, the moments may be 
calculated by using the discrete version of Eq. (2); 

X y 

However, a, region is often represented by its boundary. This Is 
the case when edge detection is used to separate objects in a scene [3]. 
It is possible to ref’onstruct the region from its boundary and use 
Eq. (5)- However, this is computationally inefficient, since a region 
usually contains many more points than its boundary. Therefore, a 
method that computes moments while traversing the boundary is desirable. 
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In this paper, such a method is derived for hoimdaries app'orntaated 
hy polygons. 


Chain code is an ordered list of numbers that represent the orient* 
tations of segments comprising the boundary. It is a ccomion]^ used 
special case of polygonal approximation vith muay attractive properties 
[U]. The method presented in this paper is used to find an algorithm 
that minimizes the time it takes to calculate moments ^om chain-encoded 
boundaries . 

In Section II, we derive the general formula for computing region 
moments from a boundary representation. The initieLl steps of the deri- 
vation rely on the theory of differential forms. Reader i. unfamiliar 
with this branch of mathematics are referred to Flanders [5]. In 
Section III, we explore the possible choices of the independent param- 
eter in the general formula to obtain equations suitable for ccm^tation. 
In Section IV, the execution time of these equations is analyzed for 
chain-encoded curves. Finally, Section V presents an algorithm that 
computes arbitrary moments from polygonal boundaries and an algorithm 
that is optimal for canputing low order moments from chain-encoded 
curves . 
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II. DERr/ATION OF THE GENERAL FORMULA 


The strategy of this derivation is to reduce the stirfaee integral 
in formula (2) to a line integral, using Stokes* theor«n. In our case, 
we want to find functions A ■ A(x,y) and B « B(x,y) such that 

m = J'x^y^dxdy » J*Adx + Bdy ( 6 ) 

R 3R 

where 3R is the boundary of R. Stokes* theorem states that for Eq. ( 6 ) 
to hold, we must have 

x^y*^dxdy * d(Adx + Bdy) (7) 

where d, the differential operator, is defined in two dimensions by 
d ■ dx^+ dy— 15]. Evaluating the d operator in Eq. (7) gives 

p q , . /3B 3A\ , , 

* y dxdy =(—.—) dxdy 


or 


( 8 ) 


p q 3B 3A 
3x 3y 

The solution to this partial differential equation is 
A = ax^y*^"*^^ , B = bx^'*^^y'^ ( 9 ) 
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where a and h eu'e real numbers. We cut find a and b by substituting the 
values of A and B ftrom Eq. (9) back into Eq. (3). This yields the 
following constraint: 

b(p + 1) - a(q + 1) - 1 (10) 

The line integral to be solved then becomes 

m ■ (aydx + bxdy)x*^y*^ (ll) 

3R 

where a and b are constrained by Eq. (10). For any piecewise continu- 
ous boundary representation, Eq. (ll) will give a region's moments 
from its boundary. 

Let 9R, the boundary of the region, be approximated by a closed 

n 

polygon as in Fig. 1. Then, 3R where the are linked oriented 

line segments with endpoints (x^^ _ Vj, _ respectively. 

Since 3R is a piecewise continuous function, the integral in equation 
(ll) can be broken into a sum of integrals over each boundary se^sent. 

Let 


m 


pql 


S (aydx + bX(J^)x^y‘^ 


( 12 ) 


Then 


m 


P<1 


n 


E 


pq£ 


(13) 
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I 

i 


as Doraerea by the dotted line, 
roximated ly an oriented polygon, 
f the links and oriented in a 


The problem is to find a by Integrating Eq. (12). Let 

pqjt 

- Xj_^^ and dy^ points (x,y) along the 

line segment can be parameterized as follows: 

X ■ dXj^t ♦ Xjj, y - dy^t ♦ y^ -1 s t s 0 (li*) 

dx ■ dx^dt, dy ■ dy^dt 

Using Eq. (l**) and the bincmial theorem. 


P <1 

r Y** B 


x*^y 


t 




i«0 


L (3) 


J-O 


(!) (?) 


(15) 


i"0 J“0 


Finally, by substituting the parameterized x, y, dx, dy, and x^y*^ of 


Eqs. (lU) and (15) into the integral that defines ®pq£ io Eq. (12), th'* 


general formula is obtained: 


>qt 


U 

/ 


t»-l 


(•♦b) AXj dy^t * ayjj AXj ♦ bx^ dy^ 


P q 


2 £ ^ y?"’’ ’■ 


1»0 J"0 




i^ere b(p+l) - a(q+l) ■ 1. Equation (l6) ii not only a foraula for the 
(p4q)-th order Booent; it also contains every possible foraula for the 
(p4q).th order Boaent. Each choice of paraiMtera a and b, subject to the 
constraint given above, generates a nev valid eqxtation for a 



III. VALUES FOR THE IHDEPEHDENT PARAMETER 


Before Eq. (16) can be evaluated, real values must be assigned to 
the parmeters a and b. The goal is to choose the a and b that will m«dte 
canputing m^^^ as simple as possible. The most obvious simplification 
is elimination of one or more of the terms in the integral. It follows 
directly from Eq. (l6) and the constraint on a and b, that only three 
parameter choices will eliminate one of the terms in Eq. (16): 

Choice 1; a ♦ b ■ 0, b ■ l/(p+q+2) 

Choice 2: a ■ 0, b ■ l/(p+l) ( 17 ) 

Choice 3; b ■ 0, a ■ -l/(q+l) 

We now look at each case in detail. For choice 1, the general 
formula reduces to 


m. 


pq£ 


btxjSyj- 


y,4x,) 


/ 

t»-l 


P 

E 

i-o j 


E(\) ( 3 ) 


■r‘ y?’’ “■‘I ‘‘A 

I ^ I I •' I 


dt 


(18) 


Noting that 




and setting the formula for choice 1 beccanes 


P q 


” *t 2 2 (') (j) 

i=0 j=0 


a. 


pqi 


(19) 


For choice 2, where a = 0 and b * l/(p+l), the middle term in the 
integral in Eq. (l6) drops out, leaving 

0 


"pqi * *’ / ‘'“n ''V * *t 


t=-l 


2- 2 (i) (j) 

i=0 J=0 


Evaluating the integral in (20) gives 


( 20 ) 


m . = b 
pq£ 


i=0 J=0 


p q 


EE";;;??)® 

i*0 J=0 


P-i+1 yq~J 


(21) 


The two double sums can be combined by making the substitution 
k=i+l in the first double siun. Making this substitution yields 
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k-»-J 


P+1 q 

k*J*l 

k»l 




P \ n\ -P->^+l „<l-J a.J^-aJ+1 


(k!i) (1) 


j/ yr 


j+i 

£ 


i"0 J"0 


( 22 ) 


Finally, using the relation 
eombinel as follows: 


P+1 q 

%«"’2 E 

i»0 j«0 


iiiiili (P«) (1) xP-‘*l „1-J A,‘ .vJ« 

i+J+1 \ i / \J/ £ 


y; " 


(23) 


Choice 3, the case where b * 0 imd a » -l/(q+l), reduces Eq. (l6) 


to 


m a * a 
pq£ 


; 

t*-i 


;4x, Ay,t ♦ y, Ax^) 


P <l 


E T . ii ) G ) 


i«0 J»0 


(CM 


Comparison of Eqs. (20) and (ClM allows us to exploit the symmetry of 
the a ■ o b ® o cases to write the solution of Eq. (2l0 as 


P q+1 


o < 


m 


pq£ 


-1 q-5+1 * i+1 A J 

y; ■’ Ay‘ 


(2M 


1«0 1»0 
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Now we have three formulas for m^^^ - Iqa. (19) » (23) » and (25) - 
corresponding to three choices of the independent parameter. Which 
formula should be used for computation? In the next section, we will 
examine the special case of chdin-encoded curves, when dx and take on 
only the values zero, one, and minus one. For now, let us consider the 
general case of polygoned approximation, where dx and dy are arbitrary 
real numbers. We assume that eill monents up to order n « p q are to 
be ccmiputed in a single boundary traversal. 

The computational requirements for m^^ are the same for all 
parameter choices. For p + q > 0, the number of terms, N, for each 
parameter choice can be read directly from the formulas: 

For choice 1, a ® -b, N = (p+l)(q+l) =pq+p + q + l 

For choice 2, a * 0, N = (p+2)(q+l) = pq + p + 2q + 2 (26) 

For choice 3, b * 0, N = (p+l)(q+2) ®pq+2p+q+2 

The equations in (26) show that choice 1 requires at least p 

additions fewer than choice 2 and q additions fewer than choice 3. 
Equation (26) also shows that choices 2 and 3 differ by p-q terms. 
Therefore, choice 2 is more efficient than choice 3 when p > q and 
choice 3 is more efficient when p < q. 

Equations (23) and (25) are homogeneous polynomials in x, y, dx, 

and Ay, of degree p+q+2. Equation (19) is homogeneous of degree p+q 

and requires one extra multiplication for the A term. Therefore, 

x» 

choice 1 requires at most the same number of multiplications needed for 
choices 2 and 3. Clearly, for the general case of polygonal approxi- 
mation to a curve, Eq. (19) is the logical formula for computing moments. 
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IV. MOMENTS FOR CHAIN-ENCODED CURVES 


Chain code is a special case of polygonal approximation. Each 
segment of the boundary connects a grid point to one of its eight near- 
est neighbors, represented by the numbers zero through seven. Fig- 
ures 2 and 3 show our conventions for labelling the chain code direc- 
tions, image coordinates, and the positive orientation of the boundary. 


Using a chain-encoded boundary representation, we wotild like to 
compute the zeroeth through second order mcanents; m^Q, ni2.o’ ”^01’ *^^20* 
m^^^, and m^g* We have a choice of two sets of formulas. The first set 
is derived from Eq. (19): 

"oo ■ 

£=1 

n 

®10 3 2 ~ 2 ^l) 

i-1 

n 

“Ol ■ 3 2 

Jl-1 

(27) 

n 

"20 ' r 2 3 '^!) 

)l=l 

n 

"11 ' vS (*£>'£ - 2 \ '"t - I i'* “t * 3 

t=l 

n 

"02 “ V Z) *1 (>'£ - '"'l * 3 '"'t ) 

t=l 
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X 1 


Figure 2. The eight possible directions bet’ 
and its nearest neighbors are rep: 
chain code numbers as shown above 



Figure 3. The region R, shown here as bordered by the dotted 
line, has its boundary, 0R, approximated by a 
chain-encoded curve. The chain code for 9R is 0,0, 
0,0,0,0,0,7,0,b,6,6,U,5,5»5,^,5,3,^,2,i*,3,2,2,2,l,l. 



The second set of fomulas viU derive from Eqs. (23) and (23) by 
using (23} if p 2 q» and Eq. (23) otherwise. The formulas are 


M 

"00 ■ 2 K ‘“'i - 5 “t “yi) 


£■1 


aa 

"10 ■ 2 £ * 3 “l 


£=1 


"01 * - 


£»1 


aa 

"20 ' 3 £ % - 2 “^1. - ¥ 


(28) 


£=1 

n 


"u " IS - I "t - ==1 


£-1 


* 3 "t “"t * 3 >^1 “‘Wt - 1 “» ‘^D 

n 

"02 ' - IS “t • 2 ='« “t 

£=1 

These formulas can be evaluated most efficiently if the terms of 

each sum are accumulated separately as the boundary is traversed. Let 

PMpq[k] represent the kth term, in the order written above, of moment 

m . For example, PM10[2] =» x.Ax.Ay, in Eq. (28). After completing the 
pq K Jo A 

traversal, each PMpq[k] term is multiplied by its coefficient. These 
quantities are then combined to obtain the moments. 
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The chain code representation allows another siaplification. Since 
and are limited to the values zero* onet and minus one, it is not 
necessary to perform the multiplications hy dx^ and dy^ indicated in the 
formulas. If dx^ terms containing dx^ ^^ 2 ^ 

ignored. Otherwise, dx^ and dy^ simply determine the sign of each term. 
The values of dx^ wd dy^^ are completely determined hy the chain code num- 
ber. We therefore treat each of the eight chain code directions as a 
separate case. Table 1 shows the eight cases for m^^, using the formula 
frcxn Eq,. (28). 


Now, we must decide which set of formulas executes in the least 
amount of time, using the implementation strategy discussed above. 

Clearly, the formulas in Eq. (28) have more terms than those in Eq. (27). 
However, this does not tell the whole story, since the restrictions 
on and dy^^^ mean that each term is not necessarily canputed every time. 
The probability that a term is computed equals the probability that its 
dXj^ or dy^^ factors are nonzero. Let be the probability that the kth 

term is computed and let P(z) represent the probability that event z 
occurs. Then, 


TT 


k 


1. 

if 

P(dxf‘0), 

if 

P(dyj^O), 

if 

P(dxdyf^O), 

if 


there are no dx or dy factors present 
a dx but not a dy factor is present 
a dy but not a dx factor is present 
both dx and dy factors are present 


(29) 
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Table 1. Chain code values and an example of moment calculations 



I 

i 


Chain code 

0 

B 

2 

3 


5 

6 

7 

dx 

1 

■ 

0 

-1 

-1 

-1 

0 

1 

dy 

0 

H 

-1 

-1 

0 

1 

1 

1 

PMOO[l]«xdy 

0 

H 

-X 

-X 

0 

X 

X 

X 

PM00[2]=dxdy 

0 

D 

0 

1 

0 

-1 

B 

1 


The eight chain code values are shown with corresponding values of 
dx, dy, and the terms PM00[l], PM00[2] used to compute m^Q. 
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Qrocn and Verbeek [6] have calculated that, for an ucbltrary 
chain-encoded figure, the probability of an even chain code Unk' la 
0.5658 and the probability of an odd link ia 0.hll»2. Therefore, refer- 
ring to Table 1, 

PCdxf^O) ■ P(odd code) ♦ §■ P(ev«a code) ■ 0.7071 

^ (30) 

P(^f*0) ■ P(odd code) + ^ P(even code) ■ O.7071 

PCdxdyftO) ■ P(odd code) ■ 0.1»ll»2 

Let be the execution time of Pt^q[k], M the execution time for 
one multiplication, and A the execution time for one additidn. Let 
be the number of multiplications needed to calcualte PMpq[k]. Then the 
expected execution time for this term, E(Tj^), is given by 

E(Tj^) ■ iT^ (y^ M + A) (31) 

Table 2 gives Uj^, and E(T^) for every tenn in Eqs. (27) and, 
(26). The total expected execution time E(T) Is Just the sum of the 
E(Tj^) over all terms in the formula. The result of adding the E(Tj^) 
in Tables 2 and 3 reveals that 

For Eq. (27), E(T) - 5M + 12.1*852A 

(32) 

For Eq. (28), E(T) - 3.9>*97M + 10.33U3A 

Therefore, ve viU use equation (28) to calculate the mcsoents from 
a chain- encoded boundary. 
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Table 2. Coaputational requiraaente for Eqaationa (27) 



k 

fMwUl 


‘'k 


®00 

1 


1 

0 

0 

“lO 

1 

V 

1 

1 

1 


2 

A,4X 

0.7071 

0 

0 

“oi 

1 

A,y 

1 

1 

1 

“20 


A-dy 

*■ 2 
A.x 
1 

0.7071 

1 

0 

1* 

0 

1 



A.xdx 
*■ 2 
At^ 

0.7071 

0.7071 

0® 
0 , 

0 

0 

“ll 


A|^xy 

1 

1 

1 



A^xAy 

0.7071 

o‘ 

0 


3 

A^ydx 

0.7071 

0*^ 

0 

CVl 

o 

B 

■ 

A AxAy 
2 

V 

O.Ull*2 

1 

0 

1 

0 

1 



A.ydy 

0.7071 

0* 

0 


H 

2 

A,dy" 

0.7071 

0 

0 

Total 



12.0710 


5 


The expected maaber of additions is ,J»l42 £n^ ■ 12.1*852.^ 

The expected nvaaber of multiplications is * 5. 


All or part of an indicated product is computed for a previous term; 
e.g., P»C0[l] ■ AjX^ ■ (Aj^x)x ■ PM10[l]»x. 

^There is one extra addition required to cosqpute A, • xdy - ydx which 
is performed with probability P(AxAy # O) ■ ,Ull*27 
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Table 3> Coaputational requirnnenta for Equation! (28) 



m 

PQ 

k 

PMpqlkl 


1 

‘‘k 

”k^k 

“oo 

1 

xAy 

0.7071 

0 

0 


2 

AxAy 

O.Ul42 

0 

0 

"lO 

1 

x^Ay 

0.7071 

1 

0.7071 


2 

xAxAy 

0.i*ll*2 

0 

0 


3 

Ax^Ay 

O.J*ll42 

0 

0 

"oi 

1 

y^Ax 

0.7071 

1 

0.7071 


2 

yAxAy 

0.i»lU2 

0 

0 


3 

AxAy^ 

0.1*lJ+2 

0 

0 

“20 

1 

x^Ay 

0.7071 

1*^ 

0.7071 


2 

x^AxAy 

0.1ilU2 

0* 

0 


3 

2 

xAx Ay 

0.1*11j2 

0 

0 


k 

Ax^Ay 

0^ 

0 

0 

“ll 

1 

2 . 
x yAy 

0.7071 

1* 

0.7071 


2 

2 2 
x'^Ay'^ 

0.7071 

0* 

0 


3 

xyAxAy 

0.1»lU2 

1 

0.141142 


k 

2 

xAxAy 

0.1al*2 

0 

0 


5 

yAx^Ay 

O.U 1 U 2 

0 

0 


6 

Ax^Ay 

O.U 1 U 2 

0 

0 

CVI 

o 

8 

1 

y^Ax 

0.7071 

1* 

0.7071 


2 

y“^AxAy 

0.i»ll42 

0^ 

0 


3 

2 

yAxAy 

0.t*ll»2 

0 

0 


1* 

AxAy^ 

0^ 

0 

0 

Total 



IO. 33 U 3 


3.9*»97 


The expected number of additions is ■ 10.33^3. 

The expected number of multiplications is ■ 3.9^97. 


All or part of an indicated product is computed f ir a previous term; 
e.g., PM20[l] » x^Ay ■ x(x^Ay) ■ x*PM10(l], 

\hese terms are not computed, since for Ax, Aye(-l,0,l},AxAy ■ Ax^Ay 
■AxAy’. Only PM00[2] ■ AxAy is computed. 
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V. THE ALGORITHie 


first algorithm prestnted hare cui^tea a^Benta from arbitrary 
polygonal boundazy curves, using the fomulas in equation (27). The 
folloving information is needed as input: 

SVERT ~ the number of vertices in the polygon. 

Xli], Y[i], i ■ 0, 1, 2, ...» NVERT — ordered lists of the x and y 
cocrdinates, respectively, of the vertices of the polygon. For 
convenience, X[0}« X[NVERT] and Y(0] ■ YtUVERT]. The coordinates 
are labelled by the convention shown in Fig 3. (0,0) is the upper 

lefthand corner of the image; the axis points to the rij^t, and 
the +y axis points downward. 

MAXORD — the highest order m<mtent to calculate 

In addition, it assvmcd that the procedure BC0EFF(m,n) which cal- 
culates the binomial coefficient has been declared. 

The output of this algorithm is a list of all (p + q)-th order 
roOTients where P ^ q MAXORD. The moments are stored in lexicographical 
order in the eurray MOMENT as follows: 

“OO’ ®10* ®01* ®20* ®11’ ''o;. ®0MAX0RD 

All entries of the array MOMENT are assumed to be initialized to zero. 

BEGIN 
INTEGER L 

FOR L:»l STEP 1 UNTIL NVERT DO 
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CC^tlBIT traverse bo\ffldury; 




IHTEOER X,Ni 

REAL m,TAX,DELTA¥^i 

DELTAX-.-XtL] - X[L-l]i 
DELTAY!»Y{L1 - YlL-11; 

AL;-X[L]»DELTAY - Y[l1*DELTAXj 
NO}IERT[l]:-IIOMBra'[l] * AL; 

K:*li CCXOIERT K indexes the array HONEST; 

FOR N:>1 STEP 1 URTZL NAXORD DO 

BEGZR CCMIERT calculate Boeents of order M; 

ZRTEGER P.Q; 

FOR P;-M STEP -1 UHTZL 0 DO 
SECn CCM)€E!ST select next p,q; 

ZinEGER Z.SZ; REAL MPQL; 

Q;^P; 

MPQL;»0; K;-K+l; SZ:-1; 

FOR Z:-0 STEP 1 UHTIL P DO 
]^GZN C(^9fENT outer sub; 

ihteger j,sj,pz 

PZ;-BCOEFF(P.Z); 

5Z:>-8Z; 8J;«SZ; 

FOR Js-0 STEP 1 UHTZL Q DO 
BEGZH CCAMEHT inner sub; 

ZHTEGER E30>: REAL T; 

C(M<ERT initialize partial product for this tens; 
T!»SJ;»-SJ; 
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FOR EXP;=I STEP -1 UNTIL 1 DO T:»T«DELTAX; 



I 


I 


FOR EXP:*J SlUP -1 UNTIL 1 DO T;=T»DELTAY; 1 

FOR EXP;=P-I STEP -1 UNTIL 1 DO T;=T*X[l]; 

FOR EXP:»Q-J STEP -1 UNTIL 1 DO T:=T»Y[L]; 

MPQL:=MPQL + (PI«BCOEFF (Q,J)»T)/l+J-l) ; i 

I 

END; COMMENT end of FOR-J loop; 

END; COMMENT end of FOR-I loop; 

MOMENTfK]:-=MOMENT[K] + AL*MPQL; 

i 

END; COMMENT end of FOR-P loop; 

END; COMMENT end of FOR-M loop; 

END; COMMENT end of boundary; 

BEGIN COMMENT scale moments by l/(p+q+2) 

INTEGER K,M,I; 

K:=0 

FOR M;=0 STEP 1 UNTIL MAXORD DO 
FOR I:=0 STEP 1 UNTIL M DO 

BEGIN K:=K+1; M0MENT[K] :=M0MENT[K] /(M+2) END; COMMENT M=P+Q; 

END; 

END; 

The second algorithm computes specific low order moments from 
chain-encoded curves. It needs the following information; 

X,y — the coordinates of the starting point of the chain. The 
coordinate labelling convention described above is used. 

L — the number of links in the chain. This is the number of 
segments that make up the boundary curve. 
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CC[i], i * 1, L - the list of chain code representing the 

curve, CC[i] e{0, .... 7) as shown in Fig. 2. The chain code 
is ordered so that the curve is traversed in a clockwise direction. 

This algorithm generates a list of low order moments stored in the 
array MOMENT as follows: 

*“00’ ®10* ®01* ®20’ “ll’ "02‘ 

BEGIN 

INTEGER ARRAY PM00[l:2],PM10,PM01,PM20,PM02[l:3] ,PMll[l:6] ; 
INTEGER I,X2,X3,Y2,Y3,XY,X2Y; 

COMMENT initialize sums; 

PM00[1]:=PM00[2]:=0; 

FOR I:=l STEP 1 UNTIL 3 DO 

PMIO [ I ] : =PM01 [ I ] : =PM20 [ I ] ; =PM1 1 [ I ] : =PM02 [ I ] : =0 ; 

PMll [ 1+ ] : =PM11 [ 5 ] : =PM11 [ 6 ] : =0 ; 

FOR I:=l STEP 1 UNTIL L DO 
CASE CC[I] OF 

BEGIN COMMENT case number corresponds to chain code value; 

[0] BEGIN COMMENT Dx=l, Dy=0; 

X:=X+1; 

Y2;=Y»Y; Y3:=Y2*Y; 

PM01[1] :=PM01[1] + Y2; 

PM02[l] :=PM02[l] + Y3 END; 



[1] BEGIN COMMENT Dx»l, Dy^l; 

X:»X+1; 

X2;»X»Xi X3:=X2*Xi Y2;«Y*Yi Y3;-Y2*Y; XY;-X*Yi 

X2Y:»X*XYi 

PM0O[l]:»B«0O[l] - X; 

PMOO[2]:*PMOO[2] - 1; 
m0[l]:=PM10[l] - X2; 

PM10[2]:»PM10[2] - X; 

PM10[3];=PM10[3] - 1; 

PMOl[l];=PMOl[l] + Y2; 

PM0l[2]:=PM0l[2] - Y; 

PMOl[3]:=P^lOlI3] + 1; 

PM20[1]:=PM20(1] - X3; 

PM20[2]:=PM20(2] - X3; 

PM20[3] :=PM20[3] - X; 

PMll(l]:=PMll[l] - X2Y; 

PM11[2]:=PM11[2] + X2; 

PM11[3]:=PM11[3] - XY; 

PMll[l*]:=PMll[i*] + X; 

PM11[5]:=PM11[5] - Y; 

PM11[6]:=PM11[6] + 1; 

PM02[1]:=PM02[1] + Y3; 

PM02[2]:=PM02[2] - Y2; 

PM02[3]:=PM02[3] + Y END; 

[2] BEGIN COMMENT Dx=0, Dy=-1; 

Y:-Y-l; 

X2;=X*X; X3:=X2*X; X2Y:=X2»Y; 
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PM00[ll:=m)0[l] - X; 

PM10[l]:»PM10[ll - X2; 

PMll[l]:*PMll[l] - X2Y; 

PM11[2]:=PM11[2] + X2; 

PM20[l]:*PM20[l] - X3 END; 

[3] BEGIN COMMENT Dx=-1, Dy=-1; 

X:=X-1; Y:=Y-1; 

X2:=X*X; X3:=X2»X; Y2:=Y»Y; Y3:=Y2*Y; XY:=X»Y; 

X2Y:=X«XY; 

PMOO[l]:=PMOO[l] - X; 

PMCX)[2]:=PM00[2] + 1; 

PM10[l]:=PM10[l] - X2; 

PM10[2]:=PM10[2] + X; 

PM10[3]:=PM10[3] - 1; 

B40l[l]:=PM0l[l] - Y2; 

PM0l[2]:=PM0l[2] + Y; 

PM0l[3];=PM0l[3] - 1; 

PM20[l]:=PM20[l] - X3; 

PM20[2]:=PM20[2] + X2; 

PM20[3]:=PM20[3] - X; 

PMll[l];=PMll[l] - X2Y; 

PM11[2]:=PM11[2] + X2; 

PMll[3]:=PMll[3] + XY; 

PM11[!+]:=PM11[4] - X; 

PM11[5]:=PM11[5] - Y; 

PM11[6]:=PM11[6] + 1; 

PM02[1]:=PM02[1] - Y3; 
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PM02[2]:=PM02[2] + Y2; 

PM02[3]:=PM02[3] - Y MD; 

[U] BEGIN COMMENT Dx=-1, Dy*0; 

X:=X-1; 

Y2:=Y*Y; Y3:=Y2*Y; 

PM0l[l]:=PM0l[l] - Y2; 

PM02[l]:=PM02[l] - Y3 END; 

[5] begin COMMENT Dx=-1, Dy=l; 

X:=X-1; Y:=Y+1; 

X2:=X*X; X3:=X2»X; Y2:=Y*Y; Y3:=Y2*Y; XY:=X*Y; 

X2Y:=X*XY 

PM00[l]:=PM00[l] + X; 

PM00[2]:=PM00[2] - 1; 

PM10[l]:=PM10[l] + X2; 

PM10[2]:=PM10[2] - X; 

PM10[3]:=PM10[3] + 1; 

M01[1]:=PM01[1] - Y2; 

PM0l[2]:=PM0l[2] - Y; 

PM0l[3]:=PM0l[3] - 1; 

PM20[l] :=PM20[l] + X3; 

PM20[2]:=PM20[2] - X2; 

PM20[3]:=PM20[3] + X; 

PMll[l]:=PMil[l] + X2Y; 

PMll[2]:=PMll[2] + X2; 

PMll[3]:=PMll[3] - XY; 

PM11[1+] :=PMll[lt] - X; 
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PM11[5]:=PM11[5] + Y; 

Bai[6]:=PMll[6] + 1; 

PM02[l]:=m02[l] - Y3; 

PM02[2]:=Bd02[2] - Y2; 

PM02[3]:=PM02[3] - Y END; 

[6] BEGIN COMMENT Dx=0, Dy=l; 

Y:=Y+1; 

X2:=X»X; X3:=X2»X; X2Y:=X2»Y; 

PMOO[l]:=PMOO[l] + X; 

PM10[1] :=PM10[l] + X2; 

PMll[l]:=PMll[l] + X2Y; 

PMll[2] :=PM11[2] + X2; 

PM20[l]:=PM20[l] + X3 END; 

[7] BEGIN COMMENT Dx=l, Dy=l; 

X:=X+1; Y;=Y+1; 

X2:=X»X; X3:=X2»X; Y2:=Y*Y; Y3:=Y2«Y; XY:=X*Y; 

XY2:=X*XY; 

PMOO[l] :=PMOO[l] + X; 

PMOO[2]:=PMOO[2] + 1; 

PM10[l]:=PM10[l] + X2; 

PM10[2] :=PM10[2] + X; 

PM10[3]:=PM10[3] + 1; 

PMOl[l] :=PMOl[l] + Y2; 

PM01[2]:=PM01[2] + Y; 

PMOl[3]:=PMOl[3] + 1; 


i 



j 
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PM20[l];»PM20[l] + X3j 
PM20[2];»PM20[2] ♦ X2; 
PM20[3]:»PM20[3] + Xj 
PM11 [i 3:®PM11[1] + X2y-, 
PMll[2]:=PMll[2] + X2; 
PM11[3]:*PM11[3] + XY; 
PM11[U];»PM11[U] + X; 
PM11[5]:=PM11[5] + Yi 
PM11[6]:»PM11[6] + li 
PM02[1]:=PM02[1] + Y3; 
H102[2l;=PM02[2] + Y2; 
PM02[3];=PM02[3] + Y 


END 


END; 

COffflCENT now combine partied svims to get the moments 

The moments are computed in the order MOO, MIO, MOl, M20, 
Mil, M02; 

M0MENT[1]:= PM00{l] - PM00[2]/2; 

M0MENT[2]:- (3»PM10tl] - 3*PM10[2] + PM10[3])/6; 

M0MENT[3]:=-(3*PM01[1] - 3«PM0l[2] + PM0l[3])/6; 

M0MENT[1*]:= (l+*PM20[l] - 6*PM20[2] + 1+*PM20[3] - PM00[2])/12; 

M0MENT[5]:* (l2*PMll[l] - 6»PMll[2] - 12»imi[3] + 8*PMll[it] 

+ U«PM11[5] - 3»PMll[63)/2!i; 

M0MENT[6]:=-(U»PM02[i] - 6*PM02[2] + 1+»PM02[3] - PM00[2])/12 
END; 


30 


rei^:rences 


1. M. Hu, "Visual Pattern Recognition by Moment Invariants," IRE 
Trans. Information Theory . IT-8, 1962, pp. 179-187. 

2. R. Wong suid E. Hall, "Scene Matching With Invariant Moments," 
Computer Graphics and Inage Processing . Vol. 8, pp. 16-2U, 1978. 

3. R. Eskenazi and J. Wilf, Low Level Processing for Real-Time Image 
Analysis , Publication 79-79, Jet Propulsion Laboratory, Pasadena, 
California, Sept. 1979. 

k. H. Freeman, "Computer Processing of Line-Drawing Images," Computing 
Surveys . Vol. 6, pp. 57-97, 197^. 

5. H. Flanders, Differential Forms with Applications to the Physical 
Sciences , Academic Press, New York, 1963 . 

6. F. Groen and P. Verbeek, "Freeman-Code Probabilities of Object 
Boundary Quantized Contours," Computer Graphics and Image 
Processing . Vol. 7, pp. 391-^02, 1978. 


NMA-JPl-CwnI . lA. CM 


31 


